How to Run Through Walls: Dynamics of Bubble and Soliton Collisions 
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A. 



It has recently been shown in high resolution numerical simulations that relativistic collisions of 
bubbles in the context of a multi-vacua potential may lead to the creation of bubbles in a new 
vacuum. In this paper, we show that scalar fields with only potential interactions behave like free 
fields during high-speed collisions; the kick received by them in a collision can be deduced simply 
by a linear superposition of the bubble wall profiles. This process is equivalent to the scattering 
of solitons in 1 -I- 1 dimensions. We deduce an expression for the field excursion (shortly after a 
collision), which is related simply to the field difference between the parent and bubble vacua, i.e. 
contrary to expectations, the excursion cannot be made arbitrarily large by raising the collision 
energy. There is however a minimum energy threshold for this excursion to be realized. We verify 
these predictions using a number of 3 -I- 1 and 1-1-1 numerical simulations. A rich phenomenology 
follows from these collision induced excursions — they provide a new mechanism for scanning the 
landscape, they might end/begin inflation, and they might constitute our very own big bang, leaving 
behind a potentially observable anisotropy. 
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I. INTRODUCTION: WHAT HAPPENS THEN? 

The recent interest in the phenomenology of the string 
landscape [l| has sparked a resurgence of interests in 
the physics of bubbles and bubble collisions. An under- 
standing of such collisions promises much: a prediction 
of observational consequences of such collisions 0t3 ; an 
insight into the interaction of non-perturbative objects 
such as domain walls in general, a more complete treat- 
ment of the bubble counting measure [3, [l^ , possibly a 
new method of scanning the landscape and many other 
wonderful things. As Coleman ll| succinctly asked in his 
seminal paper on the fate of the false vacuum "...bubble 
walls begin to collide. What happens then? Can such 
events be accommodated in the history of the early uni- 
verse?" 

The kinematics of wall collisions is well understood: 
provided one knows the nature of the incoming and out- 
going walls (if indeed walls are the only by-products), the 
wall trajectories can be obtained by energy-momentum 
conservation |13 - [l5| . Our goal, on the other hand, is to 
understand the collision dvnamics\lu^ : how does the field 
configuration evolve through a collision? 

One fruitful approach to understanding the non- 
perturbative physics involved is numerical, namely lat- 
tice simulations of bubble wall collisions. This was pi- 
oneered almost thirty years ago by Hawking, Moss and 
Stewart [17l|. Despite their relatively (by today's com- 
putational standards) low resolutiorl]], their 1-1-1 dimen- 
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sional simulations discovered an interesting result that 
was not predicted by analytic methods: the collision of 
two true vacuum bubbles creates a relatively long live 
pocket of false vacuum which eventually collapse under 
differential pressure between the spacetime regions. This 
result was not well understood or exploited until recently. 
In high-resolution, 3-1-1 dimensional, numerical simula- 
tions, it was shown that not only is this formation of 
false vacuum pocket a fairly robust effect, the right kind 
of potential can produce a new, stable bubble of a differ- 
ent vacuum (/(/], henceforth EGHL). Furthermore, such 
classical transitions can be quite elastic — coherent bub- 
ble walls form immediately after collisions with little dis- 
sipative losses in the form of scalar radiation. EGHL also 
showed that whether such transition occurs is dependent 
on both the energy of the collision, and the height of the 
potential barrier between the progenitor bubble vacuum 
and the progeny bubble vacuum. The setup is illustrated 
in Fig. [II 

In this paper, we wish to extend the work of EGHL, 
beginning with the question: what determines the size 
of the collision induced excursion? By this we mean the 
field excursion at a; = (the location of the collision), 
from immediately before to immediately after the colli- 
sion (Fig. [1} . Naively, one might think that the excursion 
depends on the collision energy — the higher the energy, 
the further the scalar field could go — as if the kinetic 
energy stored in the bubble walls can be tapped to initi- 
ate motions in field space. This kind of reasoning proves 
to be misleading, however. 

Let us try to gain some insight by simplifying the prob- 
lem. Recall that bubble walls accelerate as the bubbles 



tigate the collisions of sub-horizon bubbles during a first order 
electroweak phase transition. 
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FIG. 1. (a). A space-time diagram showing the collision of two bubbles embedded in a false vacuum sea, in a frame in which 
the two bubbles nucleate at the same time. We refer to the overlap between the bubbles as the transition region - here, the 
scalar field is generally kicked by the collision to somewhere else in field space, (b). A schematic diagram showing the potential 
for the scalar field. The arrow denotes the scalar field kick as seen by an observer following the black arrow in panel (a). This 
is the kick immediately after a collision, which is realized as long as a certain minimum condition is met. Where the field 
will roll to in the long run depends on the precise layout of the landscape in the transition region. Note that only one bubble 
vacuum is shown. The other bubble might inhabit the same or some other vacuum. Note also the field space does not have to 
be one-dimensional. 



grow in size. A high speed collision therefore also means 
a collision between two large bubbles. In this limit, we 
can treat the collision event as one between two planar 
walls, in other words it is a 1-1-1 dimensional problem. 
The collision of two bubble walls is thus equivalent to 
the collision of two kink solitons. Much is known about 
the scattering of solitons in for instance the famous sine- 
Gordon model (Appendix E)) . We will be able to deduce 
some general model independent results in the relativis- 
tic collision limit, including the case of multiple scalar 
fields. And we will find that the collision induced field 
excursion is in fact bounded, and given by a simple ex- 
pression involving the field difference between the false 
and the bubble vacua f ^fll)) . There is a minimum energy 
threshold that must be reached for this maximal field ex- 
cursion to be realized, which is derived in mill These 
analytic predictions are verified with numerical experi- 
ments in mV\ We will digress a bit in ^V] to consider 
cases where the excursion takes the field through multi- 
ple local minima. 



These collision induced field excursions open up a rich 
set of questions, including: How do they impact infla- 
tion? What are the observational signatures if we reside 
in the transition region (Fig. [ij? What happens when 
there are multiple collisions? We will explore these ques- 
tions in WIl For the most part, we focus in this paper 
on scalar field models with canonical kinetic terms and 
potential interactions, and work in flat space. Relaxing 
these assumptions will be discussed in ^^ as well. 



II. THE FREE PASSAGE APPROXIMATION - 

IMPLICATION FOR THE COLLISION-INDUCED 

FIELD EXCURSION 

We are interested in solutions to the scalar field equa- 
tion. 
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which contain bubble/domain wall configurations. Here, 
V{(j)) is a potential with two or more metastable min- 
ima. In particular, we are interested in how the do- 
main/bubble walls interact in a collision. It should be 
emphasized that while we are primarily motivated by 
cosmological applications, much of our discussion car- 
ries over to the collision of domain or solitonic walls in 
broader contexts. 

For the moment, we focus on a single scalar field 4>. 
The generahzation to multiple fields will be discussed be- 
low. We work in Minkowski space and defer a discussion 
of gravity to WII Some of our numerical simulations do 
have an expanding background, but the backreaction of 
the scalar field on the geometry is not properly taken into 
account. 

Let us follow the strategy laid out in Jl and take the 
high speed/large bubble limit, in which case the collision 
problem becomes effectively 1-1-1 dimensional: 



dU + dl^ 



dV 



(2) 



where t is time and x labels the axis of collision. 

We are intereste d in a potential of the sort schemati- 
cally shown in Fig. 2(a) The 'false' (or parent) vacuum 
is denoted by 6^ 



and the the two bubble vacua are 
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FIG. 2. (a). The top panel displays part of a landscape of local minima of the potential V of some scalar field </>. S(f> is the kick 
induced by collision. The middle panel shows the solitonic profile that interpolates between cj)^ and (f> . The bottom panel 
shows the solitonic profile that interpolates between ip^ and 4> ■ (b). The top panel depicts the scalar field </!> as a function of 
position X at some time f < in the far past, before the two solitons collide. These two solitons can be thought of as part of 
bubble walls: the bubble on the left has an interior scalar field value of (f)^; the bubble on the right has an interior scalar value 
of (f)^ ; cj> in between is the false/parent vacuum. The bottom panel shows the scalar profile soon after the solitons collide. 
What is assumed is the free passage of waveforms (Eq. (O): the top panel comes from summing / to the left of /, the bottom 
from summing / to the right of /. Note how the field value in the collision region shifts from 



t> pre-coUision to (f> + 



post-collision, with Stf) = rj}^ + (f>^ — 2(f>^. In other words, the outgoing objects are different from the incoming solitons: they 
interpolate between different field values. 



and <j>^ respectively. In other words, the collision of in- 
terest is between two walls, one spanning from (p^ to 0^, 
and the other from 0'^ to (p^ . We will treat these vacua 
as roughly degenerate, so that the wall speed is nearly 
constant. This conflicts somewhat with our picture of 
generally accelerating bubble walls, but the degeneracy 
assumption is not really crucial — much of our discussion 
carries over to the non-degenerate case with the speed u 
below taken to be the wall velocity just prior to collision. 
In the rest frame of the respective solitonic walls, the 



wall profiles are denoted by / and / (Fig. 2(a) I, i.e. they 
satisfy the equations: 
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with the respective boundary conditions /(A — > cxd) = 
0^, /(A ^ -oo) = (P^, and /(A ^ c5o) = 0^, /(A -> 
-oo) = (p^' . 

Boosting to a frame in which the two solitons collide at 
equal and opposite speed m, the wall profiles as a function 
of space and time are: 



/ 



X — ut 



Vi 



/ - 



X + ut 



^/l 



(4) 



Here, / represents the right going soliton and / repre- 
sents the left going one. When the solitons are far apart, 
say at t — >■ — oo, the solitons only interact very weakly 
with each other, and hence we can represent the total 
scalar field solution as a linear superposition of the 
two solutions / and / 
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The constant shift of —cp^ is necessary so that at 
i < (prior to collision), the scalar field (p takes the 
(false/parent vacuum) value (j)^ between the two walls, 
and the field to the far left and far right takes the values 
(p^ and (j)^ respectively (upper panel of Fig. 2(b)). 

While it is easy to see that Eq. ([5]) is a good ap- 
proximation when the two solitons are far apart, the key 
insight in understanding soliton collisions is that Eq. (0) 
remains a good approximation even during and slightly 
past the collision time t = 0, for sufficiently energetic 
collisions i.e. w — )• 1. In this limit, both the spatial and 
time gradient terms dominate over the potential term in 
Eq. ([2]) at the walls (the wall profiles effectively become 



step functions): 
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dV 



df<l> + dl^^O. (6) 



Meanwhile, away from the wahs, the scalar field sits 
roughly at the corresponding local minimum where the 
derivative of the potential is also small. 

A superposition of the two waveforms / and / as in 
Eq. ([5]) should therefore be a good solution to the linear 
wave equation ^, as long as u ^ 1. This should hold 
even post-collision. (Let us postpone for the moment a 
discussion of when this approximation Eq. [5] ceases to 
be accurate.) Taking this seriously gives a very interest- 
ing prediction: around a; = where the collision happens, 
(p = (j)^ before collision, whereas (f) — (f)^ + (f)^ —cj)'^ after 
collision. In other words, compare the upper and lower 
panels of Fig. 



2(b) 



the pre-coUision configuration is ob- 
tained by superimposing / to the left of /, whereas the 
post-collision configuration comes from superimposing / 
to the right of / (taking care to include the constant shift 
—0^ as in Eq. [5]). In other words, the two wall profiles 
(or waveforms) freely pass through each other, creating a 
region where the scalar field takes a value that is neither 
the original parent vacuum value cj) , nor the bubble inte- 



rior values 6 or 



kB' 



We say that the collision produces 



a kick to the scalar field of 
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This is illustrated in Fig. 2(b) We refer to this observa- 



tion as the free passage approximation. 

It is worth noting that even as the two profiles or wave- 
forms / and / freely pass through each other, the nature 
of the solitons has actually changed: while the incoming 
solitons interpolate between (j)'^ and (j)^ (and between 
(p^ and (f)^ ), the outgoing objects interpolate between 
(p^ + Sip and (p^ (and between 0"^ + Sep and (p^ ). Indeed, 
the outgoing objects strictly speaking are not even soli- 
tons, since p -\- 5(p is not in general a stationary point 
of the potential (see Fig. 2(a) I. One could think of the 
collision process as a scattering event. The free passage 
of waveforms is reminiscent of the impulse approximation 
used in high energy scattering. 

The free passage approximation should break down 
soon after free passage itself. As illustrated in the lower 
panel of Fig. 



2(b) 



there is a region between the out- 
going objects where the field gradient is small while the 
derivative of the potential is non-negligible (i.e. (p"^ + S<p 
is not a stationary point of the potential) . Eq. © there- 
fore ceases to be a good approximation, and one must 
account for the effect of the potential. 

Under what condition then, is the free passage kick 
Sep = (p^ + (p^ —2(p^ actually realized? In the subsequent 
sections, we will estimate, and verify with numerical cal- 
culations, a minimum collision energy /velocity in order 
for the free passage kick to be successfully executed. If 
the collision is not sufficiently energetic, the scalar field 
can still make a fleeting excursion of a similar size, but it 



quickly retreats back to the false or bubble vacua there- 
after. We will refer to this as an unsuccessful or failed 
kick. Anticipating a bit numerical experiments we will 
run, the distinction between a successful and a failed kick 
is illustrated in Fig. |6] & [T] 

It is important to emphasize that our primary goal here 
is to work out the fleld excursion shortly after a collision, 
rather than its long term behavior. As noted above, the 
free passage approximation generally breaks down soon 
after a collision (even if the kick is successful), and there- 
fore cannot be used to predict the long term outcome of 
the field dynamics. The short term outcome, however, 
does have long term implications: if the free passage kick 
fails, the field (p would inevitably roll back to the false 
or bubble vacua; if the free passage kick succeeds, the 
field's long term evolution would be driven by whatever 
the basin of attraction it happens to be in after the kick. 
If (p + 5p is in the basin of attraction of some new vac- 
uum, then we can have the formation of a new vacuum 
bubble. 

This is why the case illustrated in Fig. [5] is in a sense 
not so interesting, for no matter whether in the short 
term the free passage kick succeeds or not, in the long run 
it would not roll to any new vacuum. A more interesting 
example is one where p^ = (p^ i.e. colliding two bubbles 
inhabiting the same vacuum. This is illustrated in Fig. 
[31 

What happens in this case is that the collision induced 
kick 5<p = 2{(p^ —<p^) could lead to the formation of a new 
bubble that is at neither the parent vacuum cp nor the 
bubble vacuum (p^ . Indeed, in the example depicted in 
Fig. [21 Sep is sufficiently large that the kick takes the field 
over a new barrier (at the far left). Provided the collision 
satisfies a minimum energy condition to be derived in the 
next section, we would have an interesting outcome: the 
creation of a bubble of a new vacuum in the collision 
region. 

It is worth noting that the free passage approximation 
predicts the walls retain their integrity through a colli- 
sion (even as they change character: a wall that interpo- 
lates between (p^ and </)^ become one that interpolates 
between <p^ and cp + S(p). It is perhaps surprising that 
the collision does not result primarily in dissipation into 
scalar radiation instead. The reason is because at a high 
collision speed, the (potential) interactions are negligible 
and the field is essentially free (Eq. [6]). As we will see 
in mV[ this simple reasoning appears to be borne out 
by numerical simulations. There will be some radiative 
loss, but the amount tends to be small. The potential in- 
teractions are of course important after the free passage 
kick. For instance, in the example depicted in Fig. |31 
after a successful kick, the scalar field will roll towards 
whatever new vacuum that exists on the far left. The 
outgoing walls will adjust their profiles, as the region be- 
tween them relaxes towards this new vacuum. 

Conversely, if the collision does not satisfy the min- 
imum energy condition, then the field excursion barely 
reaches Sp = 2(0^ — (p^), and eventually retreats back 
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FIG. 3. Analog of Fig. [2] for ijf)^ = (f>^ , such that the kick S(f> — 2{(f>^ — 0'^). As before, the incoming configuration in the top 
panel of (b) corresponds to summing /(A) to the left of /(—A), while the outgoing configuration in the bottom panel of (b) 
corresponds to summing /(A) to the right of /(—A) (more precisely: the free passage solution in Eq. (O). Physically, one 
can think of this as the annihilation of a soliton and an anti-soliton giving rise to a different pair of objects. From the point 
of view of bubble collisions, we have here two (ji^-bubbles (immersed in the parent vacuum sea (j)^) colliding with each other, 



triggering a new bubble forming at 



+ 1 



(at least momentarily). 



to the false/parent or bubble vacua, (j)^ or tp^ . As we will 
see, this can happen even if, momentarily, the field has 
crossed the new potential barrier to the left. The retreat 
is due to tendency of the spatial gradient force to coun- 
teract the (free-passage implied) field excursion. In the 
language of solitons, a retreat (i.e. an unsuccessful kick) 
means the outgoing pair of (soliton-like) objects cannot 
overcome their mutual attraction, a phenomenon we will 
observe in numerical simulations (Fig. [7]). 

Lastly, to cover the range of possible outcomes, con- 
sider a situation where contrary to Fig. [3l the kick 
S(j) = 2{(j)^ — (j)^) is not large enough to take the scalar 
field over the new barrier to the left. In this case, even if 
the kick were successful, the scalar field would have to roll 
back to (p^ or even cj)^ eventually. More quantitatively, 
one can see that at the collision point (x = 0), 
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after free passage (recalling that the free passage approx- 
imation predicts vanishing space and time derivatives at 
X = after a collision). Therefore, the field retreats back 
towards 4>^ or (f>^. 

In summary, the free passage kick of Scf) = (f>^ + cj)-^ — 



2(j)^ 



or : 



) if the bubbles inhabit the same 



vacuum, constitutes a maximal field excursion shortly 
after a collision. 

Generalization to multiple fields. Let us conclude 



this section by pointing out that all of our arguments 
above translate in a straightforward manner to the case 
of multiple scalar fields, i.e. the equation of motion is 



00. ^, 



(9) 



Suppose we have a landscape of a multiple-dimensional 
scalar (j). Suppose further that we have a parent or false 
vacuum sea at cj) , within which there are two bubbles, 
one at <j>^ and one at (j>^ . When these two bubbles 
collide, the free passage approximation tells us that the 
scalar field should receive a kick in the collision region 



from 



to 



where 
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The situation is illustrated in Fig. 2] This prediction for 
5<f> should be accurate momentarily after a high speed 
collision. Note that the linear superposition (Eq. [5]) 
actually says more than this: not only is there a vector 
sum rule governing the end point of the excursion Eq. 
pO)) . the colliding wall profiles, which are not necessarily 
straight lines in field space, also add to give the excursion 
profile. 

Let us close by noting two crucial assumptions: that 
there are only potential interactions, and that the multi- 
scalar kinetic term (in the action) is diagonalized and 
canonical. We will discuss deviations from these in WII 
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FIG. 4. A schematic diagram in field space, illustrating the 
collision induced kick in a landscape of a multi-dimensional 
scalar field (f). cf denotes the false/parent vacuum value. 
Embedded inside the false vacuum sea are two bubbles, one 
inhabiting the local minimum at (j>^ , the other at cj>^ . When 
these two bubbles collide, the field in the collision region 
would undergo an excursion Sep (at least momentarily) away 
from (ji^, where S(j) is simply the vector sum of 4>^ — 4>^ and 
4>^ — (j) . The dashed lines denote the bubble wall profile 
from vacuum A to vacuum B and vacuum B' respectively. 
The solid line is a linear superposition of these two wall pro- 
files, and illustrates the field excursion upon bubble collision 
predicted by free-passage (Eq. [19|1. 



III. A MINIMUM ENERGY CONDITION FOR 
A SUCCESSFUL KICK 

Let us consider under what condition the free passage 
approximation is accurate, i.e. that the scalar field suc- 
cessfully realizes the field excursion S(f> — 2{(f>^ — (f>^) in 
the example depicted in Fig. [3] More general situations, 
such as that depicted in Fig. 01 will be considered later. 

The field configuration can be expressed as 



(j){t, x) = (t){p{t, x) + a{t, x) , 



(11) 



where 0fp is the free-passage solution which is a linear 
superposition of right-going soliton and left-going anti- 
soliton Eq. ([5]): 



0fp(i,a;) = / 
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and (7 represents the deviation from the actual solution 
from the free passage solution. The free passage approx- 
imation is accurate only when a is small compared to. 



a/5(j) < 1. We 
ut]/y/\ — V?) and 



say, the free passage kick itself i.e 
will use the shorthand fa = f{[x- 

Using Eq. Q, we can rewrite the (f> equation of motion 
as an equation for a 
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The question is: under what condition does a remain 
small after the right-going soliton and left-going anti- 
soliton have passed through each other? 
A formal solution to Eq. (fT3|) is 



a{t,x) = — 
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with the kernel function g{t' , x') 
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We estimate this integral by ignoring a in the integrand 
(first term in g) - the free-passage approximation is self- 
consistent if the resulting estimate for a is indeed small. 
We are most interested in its value at a; = 0, and by 
making a change of variables X' = x' j \/\ — v? and 
T' = ut' /yjl — u^ we can rewrite it as follows to show 
its explicit dependence on the velocity u: 



a(i,0): 
'dV 




'> = fR+fL-4>'^ 
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Note that as w — )■ 1, tr — > confirming our assertion 
that free passage is the right approximation in the rela- 
tivistic limit. The time t we are interested in is when 
the solitons have just passed through each other, i.e. 
ut/\/\ — v? = fJ.~^, where /j,^^ is the rest-frame thick- 
ness of the bubble walls. The thickness can be estimated 
by n^ ~ \d'^V/d(j>'^\ evaluated at the barrier that sepa- 
rates the parent and the bubble vacua. The integrand 
in Eq. (J16p vanishes when the two incoming solitons are 
far apart, and is dominated by when T' ^ ±/i~^. A rea- 
sonable estimate can be obtained by focusing on what 
happens post-collision, when the integrand is dominated 
by the first term dV/dcj) at = 20^ — 0"*. Therefore, we 
have 
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As we argued above, self-consistency of the free passage 
approximation demands that |cr| ^ \5(f)\ = |2(0^ — 4>^)\. 
To be concrete, we demand that for the field to continue 
marching onward after the free-passage kick, the collision 
velocity must satisfy 
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where the efficiency factor a~^ ss 5.5 is determined nu- 
merically. This efficiency factor tells us how much smaller 
\a\ has to be compared to \S(j)\ for free passage to be a 
good approximation. 

A collision that satisfies the condition Eq. (fT8)) can 
successfully realize the free passage kick. This does not 
by itself mean there is a transition to a new vacuum — 
for that to happen, the kick must take the field to within 
the basin of attraction of a new vacuum. Conversely, 
collisions that do not satisfy Eq. ([TSl) produces a fleeting 
field excursion that barely makes Sip = 2{(j)-^ — 0"*), and 
the field beats a retreat soon after. 

There exists an important caveat to Eq. (fT8|) : iidV/d(f> 
a,t (j) — 2(j)^ ~ (j)^ vanishes, this minimum energy condi- 
tion must be re-evaluated. What we have done in writing 
down Eq. ^TE\\ (and Eq. [T7]) was to approximate the 
full integrand in Eq. (|16p by one single term. This is ob- 
viously a simplification that needs to be fixed if the term 
we have chosen happens to vanish. This happens, for in- 
stance, when (p — 2<j)^ — (j)^ is another local minimum of 
the potential. 

Finally, an intriguing counter example to the above 
caveat is the well known sine-Gordon soliton, where the 
free-passage field excursion is successfully realized no 
matter what value u > takes [2l|. In this case, the 
condition Eq. (fT8|) is both trivial and exact. A brief 
summary is provided in Appendix [XI 

Generalization to multiple fields. The above rea- 
soning translates straightforwardly to the case of multiple 
scalar fields, with the equation of motion ([9]). The solu- 
tion can be expressed as = 0fp+(T, with the free-passage 
solution: 



(j)fp{t,x) = f 
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of Eqs. dH]), ([H]) and dTT]) can be written down. The 

net free-passage field excursion is given by Eq. ([TOl) . A 

crude self-consistency condition for a successful kick is 

therefore: 
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(20) 



where the derivative is evaluated at (f) = (p^ + (f)^ — 0"*, 
and /i^^ is the rest-frame thickness of the thicker wall. 



^ In this case, if the different minima are degenerate, the outgo- 
ing pair of objects are then true sohtons, and Eq. II18I I can be 
replaced by: 1/(1 — w'^) > "i^^^/m,?^, where mout and m\^ are 
the rest-mass (or tension) of the outgoing and incoming sohtons, 
i.e. m = J d(f>\/2V with the hmits of integration ranging over the 
appropriate vacua {V being set to zero at (f>^ and (p ). Energy 
conservation implies fn^^f./m?^ = {l—v'^)/(l — u'^), where v is the 
outgoing velocity and u is the incoming one. Demanding v < 1 
yiel ds the inequality. Allowing for radiation losses strengthens it 
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FIG. 5. The potential Eq. H2ip with three minima at (J)a,4>b 
and 4>c- For the 3 + 1 simulations, we nucleate two bubbles of 
pB from a sea of <J)a at the distance d = 27_Ro where Ro is the 
initial bubble radius and 7 is the Lorentz factor at collision. 



The same caveat about a vanishing potential gradient at 
the excursion point for the single field applies here as 
well. 

To close this section, let us state the transition condi- 
tion for the creation of a bubble inhabiting a new (i.e. 
neither parent nor original bubble) vacuum via a colli- 
sion: the collision energy must he high enough to satisfy 
Eq. i2U\) . and the free-passage excursion must take the 
field to within the basin of attraction of a new vacuum. 



IV. NUMERICAL SOLUTIONS 

In this section, we display a number of 1 + 1 and 3 + 
1 dimensional numerical solutions to test our analytic 
results so far, namely the free passage kick Eq. ([7]), and 
the minimum energy condition Eq. ([T5)) . We adopt the 
following potential (Fig. [5]): 



V3W = T-^- 



0o) (0 + M) +eA0g(</'-0o), (21) 



where we label the subscript 3 to signify that there are 
three minima associated with this model. There are four 
free parameters in this model, A and 4>o represent the 
overall scaling of the height and breath of the potential, 
e breaks the degeneracy of our minima and provides the 
pressure that accelerates the bubble walls, and 77 breaks 
the symmetry of the potential (making the minima un- 
evenly spaced and the barrier heights unequal). This is a 
modification of the canonical 3-minima model of EGHL 
which can be recovered by setting rj = 1. Here, we adopt 
r] — 1.15. We have neglected an additive constant to the 
potential, since we are working in the Minkowski limit. 

Let us start with results of the 1+1 dimensional simula- 
tions, for which e is set to zero. Here, we collide sohtons 
that interpolate between the appropriate minima. The 
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FIG. 6. These are snapshots of the 1 + 1 dimensional numerical solution to the coUision of a soliton and an anti-soliton, 
interpolating between 4>b and (f>A, as in Fig. (5] They are time ordered from panel (a) to (f). The red line shows the numerical 
solution. The black line shows the free passage or linear superposition prediction Eq. ^. Here, the minimum energy condition 
Eq. (|18[) is satisfied: the threshold for the potential chosen is 7crit ~ 1-7, the 7 used in the simulation is 2.3. 
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FIG. 7. Analog of Fig. [6]except that the 7 used is 1.4, and therefore the free passage kick fails, i.e. the field retreats back to 
4>A eventually. 



results are shown in Fig. [6] for a collision that meets the 
minimum energy condition. We can see that the free pas- 
sage approximation works very well, except in the very 
last snapshot. This is after the free passage kick has 
been successfully realized, the field in the collision region 
is taken to a point just shy of the new vacuum ^p. More 
precisely, the free passage kick takes the field to — 0o, 
whereas 0c is actually at — r/^Q. The field subsequently 
rolls down the potential to reach 0c — this subsequent 
evolution is not well described by free passage, as ex- 
pected. 

An analogous set of snapshots for a subcritical colli- 
sion is shown in Fig. [7l Here, the collision energy is 
not enough for the free passage kick to succeed. Note 
however that the field does make a brief excursion in the 
direction of ipc, but quickly retreats, resulting eventually 
in an outgoing pair of solitons that are just like the in- 
coming ones. The free passage approximation does not 
give the correct post collision field configuration, again 



as expected. 

For the 3 + 1 bubble (as opposed to soliton) simula- 
tions, we follow the method of EGHL for our numerics; 
the only difference being the potential is given by (PT|) 
and we neglect the expansion of the box for clarity as 
it is not essential to our work here - we have checked 
numericalljo that the results do not vary much with ex- 
pansion turned on in the presence of an overall vacuum 
energy much bigger than the energy difference between 
the minima. We find that the analytic arguments and 
1+1 soliton simulations are replicated very well in the 
3 + 1 dimensional simulations. The 3 + 1 simulations are 
crucial in confirming that there are no gross instabilities 
that might be missed in the 1 + 1 simulations. 

We nucleate two bubbles at zero wall velocity and al- 



^ We will be happy to provide the simulation results if the reader 
is interested. 
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FIG. 8. The top (bottom) panels show the time evolution of two bubbles whose centers are separated by 3.07?o (3.47?o) so that 
the bubbles achieve 7 — 1.5 (7 = 1.7) just before collision. Each diagram shows a 2D slice of the field values at equal time. 
This uses model (pTj) with 77 = 1.15 



low the pressure difference (e — 1/30) to accelerate the 
walls to the desired coUisional 7 determined by the ini- 
tial separation distance d. Figures 8(a) and |8(b)| show the 
time lapse of the results of two simulations with 7 = 1.5 
and 7 = 1.7, with a critical 7crit ~ 1.6 (slightly different 
from the value in 1 + 1 because of the difference in e). 
Comparing this to the theoretical expectation Eq. (fT8|) 
fixes a~^ w 5.5. 

In the first simulation, we can clearly see the field at- 
tempts to execute the free passage excursion, but due to 
the low collision velocity, the field does not fully transi- 
tion and retreats back into the middle vacuum (J)b- This 
final behavior is unlike the soliton case, where the field 
retreats all the way back to (pA — the difference arises 
because the soliton simulations have degenerate vacua, 
whereas 4>a has the highest vacuum energy in the bubble 
simulations. In the second simulation, the bubbles are 
nucleated at d — 3.4_Ro resulting in a 7 = 1.7 at the time 
of collision. In this case, there is sufficient kinetic energy 
in the fields during the collision to ensure a transition to 
the lowest minimum. 



V. PRODUCTION OF MULTIPLE 
WALLS/SOLITONS 

Suppose that in the course of the Scf) excursion, the 
field passes over more than one minima (and hence more 
than one barrier, see Fig. ^, then given a sufficiently en- 
ergetic collision, multiple barriers can be overcome and 
the end result is the production of multiple solitons each 
moving at a different velocity. Much of the discussion in 
the previous sections applies : the field is in free passage 
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FIG. 9. A potential with four minima at (J>a,4>b,4>c and ipo- 
Here, M, mi and m2 label the rest mass of the corresponding 
solitons. 



in the initial moments after the collision, and as the ap- 
proximation breaks down, begins to feel the potential and 
evolve accordingly. However, the field now traverses over 
more than one barrier, and if the energy of the collision 
is high enough such that the field ends up transitioning 
over more than one barrier, then the free-passage "wall" 
may split into two (or more) solitons in the aftermath. 

Nevertheless, this does not mean that these solitons 
will be long-lived: the splitting of the free passage wall 
into two or more solitons must obey both energy and mo- 
mentum conservation, and hence may result in the soliton 
with field values furthest away form the original soliton to 
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move backwards towards its mirror counterpart and an- 
nihilate each other. In the case of the sohtons sphtting 
into two, this effect can be understood by considering 
the sohtons as massive particles and the splitting as the 
decay of a massive particle (the original free passage soli- 
ton) into two less massive particles. This is then simply 
a kinetic problem which we can solve as follows. 

Consider a potential (FigE]) where the incoming soli- 
ton rest mass is M and the two outgoing sohtons have 
rest masses rai and m2- According to the free passage 
approximation, the incoming soliton passes through the 
potential barriers maintaining its field profile and veloc- 
ity. Since the field now traverses a different part of the 
potential, this profile is no longer a solitonic solution, 
thus its effective rest mass is in principle undefined. To 
make progress, we assume that this "soliton" instanta- 
neousljQ decay in the center of mass frame (denoted by 
primes) of the soliton with its incoming 4-momentum 



P^ = (M, 0,0,0). 



(22) 



Meanwhile, the two outgoing sohtons possess the follow- 
ing 4-momenta 

Pf = (7imi,7iTOiUi,0,0) , P!^ = (72TO2,72m2M2, 0,0). 

(23) 
Conservation of energy and momentum P'' = P!^ + P^" 
then allows us to solve for u'l and u'2 , which we can then 
transform back into the center of collision framqfl. As- 
suming then the first soliton mi is formed with velocity 
moving away (defined to be positive) from the center of 
collision frame, then the second soliton has the following 
velocity 



U2 



1 -f UnU 



(24) 



which must be > for the second soliton-antisoliton pair 
(mass 7712) not to self-annihilate. 

To be specific, consider the following modification to 
our toy potential Eqn. (PTjl 



^4(0) = 1^3(0) -^exp 



'-M' 



(25) 



where (j)a defines the location of an additional metastable 
minimum, and 5^ defines the width of that minimum. 



** This assumption is in keeping with the spirit of the free pas- 
sage approximation, although one can argue that the decay is 
a 2-step process in the following sense — the incoming solitons 
pass through each other in free passage, and then either gain or 
lose mass due to the fact that their profiles no longer traverse 
the original potential resulting in a change in velocities, before 
decaying. 

^ In this example, we have assumed that both incoming solitons 
have identical rest masses, hence this frame (where the incoming 
velocities are equal and opposite) is also the center of mass frame 
for the total system. 



Although 5 looks like a free parameter in this model, we 
fix it so that when e = our minima are all degenerate, 
hence, 



5 = ^K 



0o) (0a - (jiori) . 



(26) 



This potential is equivalent to the soliton potential 
shown in Fig. [HI For a model with 77 = 1.15, (pa — O.650o 
and b^ — O.80Q, one can calculate that for U2 = (i.e. the 
critical splitting velocity) 7 = 2.3. Colliding the solitons 
at 7 = 2.6, we can use Eq. ([24| to find that the second 
soliton will have a velocity of U2 ~ 0.2 and hence a split- 
ting will occur, a result which is numerically confirmed 
in Fig. [TUl 

As a final numerical test, we show that such split- 
tings occur even in 3 -I- 1 dimensional simulations. We 
use the same potential (|25p as in the solitonic case, 
except that we have added in a small linear tilt term 
(1/30) A(0 — 0^)00. Since we anticipate that a highly rel- 
ativistic collision is necessary to produce multiple domain 
walls, we begin by nucleating two bubbles at a separa- 
tion of d = 5.2Po- Figure [TT] shows the time evolution of 
these events. We can clearly see the free passage of the 
field during the collision, and then the creation and ac- 
celeration of two domain walls, making a double bubble. 

Finally, we comment on an interesting possibility. Con- 
sider the case such that the collision 7 is barely insuffi- 
cient to satisfy Eq. (|24| . In this case, a fully formed 
soliton-antisoliton pair between (pc and (po are formed 
but self-annihilate after at while. This secondary colli- 
sion can be treated as any other collision - if the right 
conditions are satisfied then we can form a pocket of (pB 
vacuum, creating a "bubble within a bubble" scenario. 



VI. DISCUSSION 

To briefly summarize, we find that linear superposition 
of wall profiles provides a good approximation to what 
happens at a high speed collision. The free passage of 
wall profiles implies a field excursion in the collision re- 
gion that is given by Eq. (ITU)) . For this free passage 
kick to be successfully realized, the collision must exceed 
a certain minimum energy threshold given by Eq. ()18|) . 
Beyond this threshold, the size of the kick itself is in- 
dependent of the collision speed, i.e. one cannot obtain 
arbitrarily large field excursions by pumping up the col- 
lision energy. When the minimum energy condition is 
met, the free passage kick can cause a transition into a 
new vacuum, if the kick takes the scalar field to within 
the new vacuum's basin of attraction. Interesting possi- 
bilities arise if the excursion traverses over several new 
vacua, and they can be understood by similar arguments 
( W\ . We have verified these statements using numeri- 
cal computations of soliton collisions in 1 -I- 1 and bubble 
collisions in 3 + 1. 

These findings have many interesting theoretical and 
observational implications. Let us discuss some of them. 
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(a) 



(b) 



FIG. 10. The interaction of a soliton-antisoliton pair in the presence of the potential denoted in Fig. [9] and Eq. (|25p with 
T] — 1.15, 4>a = 0.65(j>o, b^ — 0.8(/f>o- The incoming sohton-antisoliton pair spans the barrier from (f>A to 4>b (left figure). Given a 
sufficiently energetic collision - 7 > 2.3 for this potential - the field can transition to (po, forming a pair of solitons (and a pair 
of anti-solitons) in the process (right figure). In general, the outgoing pair of solitons possess different velocities depending on 
their rest masses and the actual collision energy. For some potentials, the process can be thought of as annihilation of the two 
incoming solitons to form a pair of temporary free passage solitons with the same mass, and then the decay of these solitons 
into two pairs of less massive solitons. Since the interaction is not completely elastic, scalar radiation visible as superimposed 
small perturbations is emitted. 
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FIG. 11. The time evolution of two bubbles whose centers are separated by 5.2Ro so that the bubbles achieve 7 = 2.6 just 
before collision. This uses the model (|25p with ri = 1.15, 4>a ~ 0.65</>o, b^ = 0.84>o- Periodic boundary conditions produce 
multiple images of this collision in these panels. 



A. Including Gravity 

Our simulations used a flat, non-expanding back- 
ground. To study a classical transition in the cosmo- 
logical context, first we need to know the necessary mod- 
ifications once we include gravitational effects. 

Since the collision happens on a short time scale, the 
corrections from gravity will not be in the collision itself, 
but in the domain wall motions before and after the col- 
lision. Without gravity, it appears that arbitrarily high 
barriers can always be traversed via collisions, as long 
as the incoming boost 7 is large enough. In an expand- 
ing background, however, there is a maximum boost for 
every collision. 



7Ma 



[RoHaY 



(27) 



where Rq is the initial size of the bubble B, and Ha 
is the bubble expansion rate for the parent vacuum A. 
Therefore, it is easy to generalize our results to an ex- 
panding background, as long as we maintain an addi- 
tional condition : the critical 7 to make the transition in 



flat space has to be smaller than 7Max in the expanding 
background, otherwise there will be no transition. 

After the transition, there are more complications. 
For example, if Vc > Vb, the domain walls can turn 
around and may or may not collide with each other again, 
thereby sealing the vacuum C region. Also when the do- 
main wall is heavy, niBC < \Vc ~ Vb\, gravity allows it 
to accelerate and eventually run away from both sides. 
These behaviors are hard to keep track of in lattice sim- 
ulations. An analytical study on these questions can be 
found in m. 



B. Cosmology After a Collision 

Many recent studies discussed the possibility of ob- 
serving signals from bubble collisions in our past. Sev- 
eral cases are summarized in [^. One important issue 
is whether slowroU inflation can be maintained after a 
violent disturbance like a bubble collision. The scalar 
field behavior we study here provides useful intuition to 
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address such a problem. 

If the field range of slowroll on the potential is super- 
planckian, which is usually known as a large-field slowroll 
inflation, then it is quite stable against disturbances as we 
do not expect the vacuum separation (which determines 
the collision induced kick) to be super-planckian. A suf- 
ficiently large flat region in the potential post-collision 
would provide a buffer for slowroll to continue unscathed. 
On the other hand, if the slowroll range is smaller, then 
it is more delicate. One example in [2^ seems to suggest 
that a collision always ends small-field slowroll inflation. 
Here we would like to argue the opposite — Having small- 
field slowroll inflation after a collision is as plausible as 
having it after a Coleman- De-Luccia tunneling. 

We start from Figure 7 in 2], which the authors in- 
terpreted as the collision pushing the field from slowroll 
directly to the end of inflation. The free passage tells 
us that right after the domain walls cross each other, the 
field ends up within a region (possibly very) near the par- 
ent vacuum (see Figure 7 of [2] ) . Because of the pressure 
difference, such a region undergoes oscillations and re- 
collisions as in 17[. In the figure cited above, the first 
two collisions are visible but the rest of them are too 
small to be resolved, forming an effective domain waljj. 
The pressure difference between the bubbles on left and 
right accelerates this effective domain wall to the right. 

The creation of this domain wall is highly dissipative 
- scalar radiation is emitted from the point of initial col- 
lision and subsequent re-collisions, and this scalar radi- 
ation is large enough in amplitude to disrupt the small 
field inflation, causing the inflaton to topple off its deli- 
cately balanced potential and thermalize. In other words, 
the failure of small-field slowroll inflation appears to be 
a result of the disturbances of the repeated re-collisions 
hidden in the effective domain wall. 

Nevertheless, as we have shown in this paper, counter 
to the intuition that collisions are violent, classical tran- 
sitions can be exceptionally gentle (recall the fairly ho- 
mogeneous post-collision region in our 1 -f 1 and 3 + 1 
simulations). This gentleness might live in harmony with 
the delicate small-fleld potential. Consider a potential in 
Fig. [121 If we arrange that (pA — (pB = <Pb — 't'Ci then by 
free passage, a collision between two bubbles of B nat- 
urally starts slowroll inflation in 4>C- One might argue 
that such a setup is finely tuned. Let us ask a differ- 
ent question: Can a Coleman-De-Luccia tunneling from 
vacuum B start slowroll inflation in 0c? The answer is 
yes, if we arrange the potential correctly - it is also flnely 
tuned! 

Although there seems to be an additional parameter 7 
involved in collisions, we remind ourselves the following 
two facts. 

• Our simulation shows that once the transition is 
allowed, increasing 7 further has no impact on the 
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FIG. 12. A potential with vacua A, B and a small field slowroll 
potential around C. 



held value right after the transition. Namely, the 
initial condition for starting the slowroll is not sen- 
sitive to the incoming boost. 

• In a multiverse we have an inflnite number of col- 
lisions with different 7's, so it is naturally scanned 
and does not require further flne-tuning. 

Therefore, all the flne-tuning is on the potential itself. 
We conclude that starting slowroll through collisions is 
not more flne-tuned than starting it through a single tun- 
neling Q 

Let us now turn to the question of the spacetime geom- 
etry inside a coUisionally formed bubble. The motivation 
is that we might perhaps live in it, i.e. a collision event 
might have constituted our beginning, our big bang. As- 
suming that there is slowroll inflation inside the bub- 
ble, we can approximately describe the geometry as pure 
de Sitter, whose general metric is given in [l3| : 

ds^ = -— - + fiz)dx^ + z^dHl , (28) 



where 



f{z) = l-- + H'z 



(29) 



When the domain walls are highly boosted, the con- 
servation of stress-tensor gives uCl 



fA[Zc)fc{Zc) = [fsiZct 



(30) 



where Zc is the collision radius of the two-dimensional 
hyperbolic section if 2- The A,B,C labels follow the la- 
beling of vacua in the rest of this paper. Neglecting radi- 
ation loss (which seems to be a good approximation from 



We have reproduced this result in both 3-1-1 and 1+1. 



'^ These fine-tunings are only necessary in small field infiation mod- 
els. In large field models it will be equally generic to start infla- 
tion through either tunneling or collision. 

* This is exactly true for massless(null) domain walls [l3l . and 
approximately true for light or highly boosted domain walls. 
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our simulations) , M = Q in both /a and Jb ■ As arranged 
in the potential Va> Vb > Vc, M in fc is also small and 
the M/z term soon becomes unimportant. The metric in 
region C is thus approximately 



ds^ = — 



dz^ 



1 



Hlz^ 



(1 + Hl:z^)dx 



z^dH^ 



(31) 



When there is slowroll inflation in region C, this will 
not be the proper slicing to describe the universe. Re- 
gion C is the forward lightcone of an H2 with radius Zc- 
Following the usual prescription — in the case of a single 
bubble — for deriving an open de Sitter slicing cover- 
ing the future lightcone of a point in the global de Sitter 
space, we can similarly derive the "open" slicing of the 
universe in region C: 



ds^ = 



^de + '-^^^^J^de 



H^ 
^c 



(32) 



Zc cosh Hct 



smh He t 
He 



I + H^jz^coshn dHi 



We can see that the metric preserves only the symme- 
try of H2, and is both anisotropic and inhomogeneous. 
The interesting trait is that both its anisotropy and in- 
homogeneity are correlated with its negative spatial cur- 
vature. Since the universe is roughly isotropic and ho- 
mogenous to 1 part to 10^ |23i|, if we are to live in such a 
bubble, there must exist a sufficient amount of inflation 
in this new universe after its birth through a collision 
[24| . While the possibility of statistical anisotropy have 
been discussed, at least at the level of the initial per- 
turbations |25j . our model predicts an interesting novel 
anisotropy in the form of an anisotropic spatial curvature, 
i.e. there exist two different curvature scales in different 
directions, an intriguing possibility we plan to explore in 
future work. Note that an anisotropic curvature shifts 
the Doppler peaks in the microwave background on all 
angular scales, making the effect easier to observe than 
most other large scale anomalies, as recently emphasized 
by [2a, I23I in the context of other models. 



C. Additional Ingredients to the Field Dynamics 

There are several possible complications to the field 
dynamics beyond what we have discussed. First of all, 
there can be derivative interactions. Our analysis should 
hold as long as the collision process does not probe energy 
scales above the cut-off for these derivative interactions, 
say A. This means the boosted thickness of our bubble 
walls 1/(7/^), where 1/ fi is the rest-frame thickness, must 
be larger than 1/A. This implies an upper limit to the 
collision speed that we can consider 7 < A//i. When 
derivative interactions are important, linear superposi- 
tion such as we have used in Eq. ([5]) no longer works. 

Another possibility is that the field metric (in the 
multi-field case) is non-trivial, and therefore the excur- 
sion trajectory is modified. A third possibility is that 



some other field is coupled to our bubble scalar field, and 
the coupling is such that this other field becomes mass- 
less in the course of an excursion. Massless particles get 
producecjj, and could significantly modify the excursion 
trajectory. We plan to investigate these intriguing possi- 
bilities in the future. 



D. Scanning of the Landscape 

The collision induced excursion offers a new mecha- 
nism for scanning a landscape of many vacua. Our results 
here touch on two aspects of this scanning. One is that 
even very high barriers can be overcome by a collision as 
long as the collision is relativistic enough (but subject to 
constraints from expansion, and from derivative interac- 
tions). It could well be that certain vacua surrounded 
by high barriers are more likely to be populated by col- 
lisions rather than direct tunneling. A second important 
feature of this scanning is that it is bounded — the colli- 
sion induced excursion does not become arbitrarily large 
by raising the collision speed. Rather, it is controlled by 
a simple vector sum rule (Fig. S]) , which tells us that the 
excursion can only be as large as the field difference be- 
tween the parent and bubble vacua. In theories in which 
tunneling between vastly separated vacua is common [29| , 
large collision excursions are also possible. Finally, this 
excursion has a specific direction, and it is interesting to 
ask what is the typical basin of attract there. In[30| is 
it suggested that due to the universal dilatonic runaway 
direction in string theory inspired models, classical tran- 
sition can lead to decompactification of extra dimensions. 



E. Multiple Collisions 

The free passage approximation and the resulting vec- 
tor sum rule in section |lT] is a simple fact of the field 
dynamics. It applies to any solitonic objects interpolat- 
ing between local minima. We can arrange a potential in 
which a collision is followed by further collisions, leading 
to a nested set of classical transitions. A spacetime struc- 
ture with multiple, overlapping collisions follows simply 
from the vector sum rule, as shown in Figure [T51 

A particularly intriguing version of this occurs when a 
bubble is nucleated in the presence of a compact dimen- 
sion |3ll ]. The bubble grows, wraps around the compact 
dimension, and eventually collides with itself. A new vac- 
uum opens up between the outgoing pair of walls. These 
walls eventually collide after traveling through the com- 
pact dimension. Further transitions and collisions follow, 
as far as the potential allows. This provides a novel way 
to realize an old idea by Abbott [32], by classical transi- 
tions rather than tunneling. 



® In[28[ it is treated as a direct generalization of preheating, but 
we remain conservative about their results. 
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FIG. 13. A spacelike slice through multiple, overlapping col- 
lisions. Following the vector sum rule, we have: 

CJ>C — 'i4>B — 4>A, <t>D = 3(j}B — ^(jyA- 



In the context of a scalar field, such a cascade of clas- 
sical transitions seems to rely on the existence of roughly 
evenly spaced vacua in the potential landscape. Interest- 
ingly, in [33| it is shown that in a model with multiple 
vacua constructed from extra dimensions and fluxes, a 
classical transition is the most natural result, as shown 
in Figure[T4l In these models, the multiple classical tran- 
sition structure in Figure [13] should be taken seriously. A 
similar but simpler structure was studied in j34| and sug- 
gested the possibility of a CFT description for an eternal 



inflating spacetime. 



ACKNOWLEDGMENTS 

We would like to thank Richard Easther, Brian Greene, 
Matt Johnson, Tommy Levi, and especially Matt Kleban, 
Alberto Nicolis and Erick Weinberg for many useful dis- 
cussions. This work is supported in part by the DOE 
DE-FG02-92-ER40699, and by the NASA Astrophysics 
Theory Program 09-ATP09-0049. LH thanks members 
of the CCPP at NYU and the IAS at Princeton for 
their fabulous hospitality. EAL thanks the hospitality 
of Kenyon College where some of this work was done and 
acknowledge the support of a Foundational Question In- 
stitute Mini-Grant. Research at the Perimeter Institute 
for Theoretical Physics is supported by the Government 
of Canada through Industry Canada and by the Province 
of Ontario through the Ministry of Research & Innova- 
tion. 

Appendix A: The Sine-Gordon Example 

The sine-Gordon potential is V^ = 1 — coscj) such that 
the equation of motion is —dfip + dlcj) = sincf). The exact 
one-soliton solution is 



4 tan ^ exp 



X — ut 



(Al) 



This interpolates between = 27r to the far right and 
(/i = to the far left. The exact soliton-anti-soliton pair 
solution is 



4 tan 



1 sinh(— uf/\/l — u^) 
u cosh(a;/Vl - u^) 



(A2) 



This solution describes the collision of an incoming (0/27r, 
2tt/0) pair, resulting in an outgoing (0/ — 2tt, — 27r/0) 
pair, irrespective of the size of u. Note that w = is not 
a solution ~ there exist no stable static solution of more 
than one sine-Gordon soliton. 
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